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Abstract 

The prevailing view in the gophysics community is that the electrical conductivity structure of the 
Earth’s continental crust over the 5-35 km depth range can best be understood by assuming the 
presence of intergranular fluids and/or of intergranular carbon films. Based on single crystal studies of 
melt-grown MgO, magma-derived sanidine and anorthosite feldspars and upper mantle olivine we 
present evidence for the presence of electronic charge carriers, which derive from peroxy defects that 
are introduced during cooling, under non-equilibrium conditions, through a redox conversion of pairs 
of solute hydroxyl arising from dissolution of H 2 O. The peroxy defects become thermally activated in 
a 2-step process, leading to the release of defect electrons in the oxygen anion sublattice. Known as 
positive holes and symbolized by h’, these electronic charge carriers are highly mobile. Chemically 
equivalent to O in a matrix of O 2 they are highly oxidizing. Being metastable they can exist in the 
matrix of minerals, which crystallized in highly reduced environments. The h’ are highly mobile. They 
appear to control the electrical conductivity of crustal rocks in much of the 5-35 km depth range. 
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42 I Introduction 

43 For many geophysical processes in the Earth’s crust it is crucially important to understand the nature 

44 of the charge carriers, which control the electrical properties of rocks. In this report we focus on 

45 conditions in the Earth’s continental crust over the 5-35 km depth range, where most major crustal 

46 earthquakes occur. Our goal is to identify the electric charge carriers that are present or can be 

47 activated in this depth range. 

48 Any cursory inspection of geological maps reveals how heterogeneous the crust can be. Magneto- 

49 telluric (MT) soundings have provided ample evidence for the electrical heterogeneity as a function of 

50 depth. However, they have also identified high conductivity zones that seem to consistently exist at 

5 1 certain crustal depth levels, as well as large-scale inhomogeneities in the upper mantle. 

52 The two most widely quoted causes for variations in the electrical properties of rocks are 

53 interconnected networks of aqueous fluids and/or of intergranular carbon or graphite films [Glover 

54 1996], In this report we point to the presence of electronic charge carriers in minerals and rocks, which 

55 have not been fully recognized in the past. 

56 Any electrical charge transport in rocks can be attributed to one or more of the following three 

57 contributing causes: (i) electrolytic conduction, (ii) electronic conduction, or (iii) ionic conduction. 

58 Electrolytic conduction will occur in rocks where aqueous fluids fill pores and the intergranular space, 

59 forming an interconnected network. Ionic transport is negligible at the temperatures most relevant for 

60 the crustal depth range of interest here but surely becomes important at high temperatures, both in the 

61 solid state and with the on-set of partial melting. Electronic conduction can occur over the entire 

62 temperature range of interest here, either through intergranular networks or through the bulk. 

63 

64 1.1 Emphasis on Intergranular Fluids and on Carbonaceous Films 

65 Aqueous fluids certainly exist in the shallow crust in the form of saline solutions filling the omni- 

66 present pore space in rocks. Open pores can form a communicating network through which fluids are 

67 able to flow under the influence of tectonic pressure gradients. Such forced fluid flow through narrow 

68 pores will produce electrical signals, known as streaming potentials. They arise when cations are 

69 preferentially retained along the walls of narrow conduits while anions are swept along with the fluids 

70 [. Dobrovolsky et al., 1989], Streaming potentials are the foundation of the electrokinetic effect, widely 

7 1 used in prospecting of the shallow crust by magneto-telluric and other techniques as well in attempts 

72 to understand pre-earthquake electrical or electromagnetic signals \Merzer and Klemperer, 1997; 

73 Mitzutani et al., 1976; Takahashi et al., 2007; Teisseyre, 1983; Yoshida et al., 1998], 

74 Deeper into the crust, as the overload of the rock column increases, pores that communicate with the 

75 surface must close as the difference between lithostatic and hydrostatic pressures becomes too large. 

76 Closing occurs around 3-5 km. Below this cap, fluids are at lithostatic pressures. They can either be 

77 trapped in situ in sediment-derived formations, which had been transported downward through 

78 tectonic processes such as faulting [Yardley et al., 2000], or be supplied from below through degassing 

79 of magmas during crystallization or through mineral dehydration above a mantle plume [Vanyan and 

80 Gliko, 1999] or above a subducting slab [ Hyndman , 1988], The presence of an electrolytically 

81 conductive intergranular phase should have a profound effect on the electrical conductivity structure. 

82 Two views emerged in the literature regarding the nature of intergranular films. One view emphasizes 

83 the presence of “water”. The other view emphasizes the presence of carbon films precipitated in the 
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84 intergranular space out of C0/C02-bearing fluids within the stability field of graphite [ Constable and 

85 Duba , 1990; Nover et al., 2005; Shankland et al., 1997]. 

86 The inferred presence of water has been widely invoked as an explanation for increased electrical 

87 conductivities in the crust and the upper mantle [ Frost and Bucher, 1994; Fuji-ta et al., 2004; Iio et 

88 al., 2002; Karato, 201 1; Khan and Shankland, 2011b; Nekut et al., 1977; Qin et al., 1992; 

89 Wannamaker, 1994] or to explain high conductivity layers in the crust at a certain depth ranges, 

90 around 20 km [Chen and Chen, 2002; Fuji-ta et al., 2004; Lee et al., 1983; Patro et al., 2005; Qin et 

91 al., 1992], The zones of high electrical conductivity seem to correlate with zones of low velocity for 

92 seismic waves, though this correlation has been called into question [ Bai et al., 2002], The presence 

93 of water has also been invoked to account for high electrical conductivity along seismically active 

94 faults such as the Parkfield section of the San Andreas Fault in California [ Unsworth et al., 1999], 

95 Intergranular carbon, carbonaceous or graphite films have been considered as an alternative to water 

96 [Nover et al., 2005; Olhoeft, 1981; Shankland et al., 1997], Carbon films have been invoked in the 

97 context of zones of high crustal conductivity around 20 km, the presence of which suggests conduction 

98 mechanisms that are not activated in at shallower depth and disappear at greater depth [Glover, 1996; 

99 Haak and Hutton, 1986; Hyndman et al., 1993], 

100 Surface analytical techniques such as x-ray photoelectron spectroscopy (XPS), secondary ion mass 

101 spectrometry (SIMS) or ion probes have indeed provided evidence for carbonaceous surface films a 

102 few atomic layers thick after rocks had been fractured in 95% C0 2 +5% CO or 5% CH 4 gas mixtures 

103 at temperatures between 350-500°C [Roberts et al., 1999], Such films are assumed to have been 

104 deposited in situ from fluids containing reduced carbon species such as CO or CH 4 that would have 

105 percolated through the rocks. Experiments with thin carbon films along the contact between two 

106 quartz crystals at 1 GPa have demonstrated that, when the temperature was raised to about 630°C, the 

107 graphite films disintegrate and disappear [Yoshino and Noritake, 2011], 

108 In the petrological literature a strong case has been made that the stable lower crust should be “dry” 

109 and devoid of percolating or stagnant fluid phases [Yardley and Valley, 1997], The argument is that, 

1 10 given the prevailing high temperatures, any “water” in these deep environments will be consumed 

1 1 1 rapidly - on a geological time scale - by metamorphic reactions leading to hydroxyl-bearing minerals 

1 12 such as amphiboles and mica, both of which are characterized by low electrical conductivities 

113 [Litovchenko and Mazykin, 1984; Meunier et al., 1983; Wang et al., 2012], The case for a “dry” lower 

114 crust has sparked an insightful discussion about alternative views [Wannamaker, 1999; Yardley and 

115 Valley, 2000], 

116 

117 1.2 AC versus DC 

1 1 8 Numerous laboratory studies have been conducted to characterize the electrical properties of minerals 

119 and rocks. Most often the electrical conductivity was measured using ac methods. The reason for 

120 preferring is that ac measurements can (i) provide more parameters than dc techniques and (ii) are 

121 better suited for experiments under the high temperatures and high pressures needed to explore the 

122 electrical properties of minerals and rocks under lower crust and upper mantle conditions [Fuji-ta et 

123 al., 2004; Liu et al., 1999; Nover, 2005; Wang et al., 2012; Yang, 2011], 

124 When an ac electric field is applied, a marked frequency dependence is often observed, which is 

125 caused by the complex interplay between conduction and displacement currents. Displacement 

126 currents from the polarization of electrons, bonds, mobile ions, defects that either have a dipole or 

127 develop a dipole in the alternating electric field, charges accumulated along grain boundaries or 
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charges trapped and detrapped at the interfaces between the sample and the electrodes. 

To appreciate the complexity of the situation we look at the resistivity p [Q cm] or its reciprocal, the 
conductivity a [Q _1 cm' 1 ], and at the dielectric permittivity 8 [F cm' 1 ]. For a linear system we can 
write the conduction current density J to tsi at time t and the electric displacement D caused by an electric 
field E: 

Jtotai = ctE + ehE/ht [1] 

and 

D = 8E [2] 

The relative dielectric constant 8 r of the medium is defined as 

e r = 8/8 0 [3] 

where 8o = 8.854 x 10' [F/m] is the dielectric permittivity of free space. Conduction currents and 
displacement currents will be out of phase with each other. 

The parameters that are actually measured are not the true values but the effective values given by 

e eff (oo) = 8 ’(oo) 008 ” [4] 

and 

S e ff (oo) = 8’(oo) - i(oo)<3”(oo)/oo [5] 


The ratio of the imaginary part to its real part defines the phase angle 'F given as: 

tan 5 = [(o8”(oo)]/a’(co) [6] 

(a) (b) 




1.5 2.0 2.5 3.0 3.5 

1000/K 


Figure la: Resistivity p, loss tangent, tan 8, and 
dielectric constant £ as a function of frequency 
co. The “true” conductivity o defined for co = 0. 


Figure lb: Example of the conductivity a of a 
pyroxene as a function of temperature from 25 to 
400°C and 0 to 10 6 Hz [Olhoeft, 1976]. 


Figure la illustrates how the resistivity p, the loss angle tan 6, and the effective dielectric constant s e ff 
depend upon the frequency oo. The “true” conductivity a is defined only at 0 Hz. Figure lb shows the 


4 





149 conductivity a of a pyroxene between 25°C to 400°C at frequencies oo between 0 Hz and 10 6 Hz 

150 \Olhoeft, 1976]. A frequency-dependent loss feature is seen, which shifts to higher temperatures with 

151 increasing frequencies but also becomes increasingly indistinct. 

152 Data as depicted in Figure lb have diagnostic value for the interpretation of laboratory experiments. 

153 However, they are not directly useful in the analysis of many sets of field data collected by MT 

154 techniques. The reason is that, while the first few hundred meters below the surface with EM waves in 

155 the 1-5 kHz range, low to ultralow frequencies in the Hz, mHz and pHz ranges are needed to penetrate 

156 deep enough into the Earth to send back information about the electrical structure of the deeper 

157 regions of the Earth [Cagniard, 1953], 

158 In order to build the case for dc measurements, it is opportune to take a step back and look at the 

1 59 electronic conductivity of materials from a broader perspective. 

160 

161 Electronic Charge Carriers 

162 Materials can be broadly classified into three groups: metals, semiconductors, and insulators. The 

163 three groups differ in the position of their valence bands (VB) and conduction bands (CB) relative to 

164 each other on the energy scale. 

165 In metals, VB overlaps with CB allowing electrons from the VB to flow over into the CB. The CB 

166 electrons are not bound to any particular atom, meaning that their wavefunctions are delocalized. 

167 Semiconductors have a modestly wide band gap, E g on the order of 1 eV. Such gaps allow some 

168 electrons from the VB to become thermally activated to the CB, where they act as negative charge 

169 carriers e’. Around 300 K the thermal energy kT is ~25 meV, where k is the Boltzmann constant, 

170 8.618 x 10 5 eV/deg, and T the absolute temperature. Assuming a Boltzmann distribution, the 

171 probability for VB electrons to be thermally activated across E g is e _Eg/kT . For a 1 eV band gap this 

172 means e' 40 ~ 4 x 10" 18 at 300 K. Electrons that make the transition to the CB leave behind defect 

173 electrons in the VB. The defect electrons are called holes, designated by h\ They carry a +1 charge. 

174 Holes in the VB are less mobile than electrons on the CB. 

175 In insulators, E g is wide to very wide, often 5 eV or more. Thus the probability for electrons to be 

176 thermally activated from VB to CB becomes vanishingly small. Numerically, for e _Eg/kT ~ 5 eV, the 

177 probability at 300 K or room temperature is on the order of 10~ 87 and still as low as 10" 29 for 900 K or 

178 ~630°C. Thus, the electronic conductivity of insulators should be zero. Any observed conductivity 

179 has to arise from the presence of defects and/or impurities. 

180 Eq. [7] links conductivity to the number density of charge-carrying species n and their mobility p: 

181 o = npzV/(kT) [7] 

182 where z is the valence of the conducting species, e the unit charge, k the Boltzmann constant, and T 

183 the absolute temperature. 

184 If the conductivity is thermally activated with an activation energy E ac t, we have: 

185 o = a 0 exp [-E act /kT] [8] 

186 where the pre-exponential factor a 0 is the product of several contributions including p, z and e. If 

187 more than one process control the conductivity, each with an activation energy Ei, E 2 . . ., we have: 

188 o = a, exp [-Ei/kT] + a 2 exp [-E 2 /kT] + . . . [9] 

189 Plotting log a versus 1/T eq. [8] will give a straight line, called an Arrhenius plot, with the slope E act . 

190 Eq. [9] will give two or more straight sections corresponding to Ei, E 2 . . . 

191 In semiconductors electrons in the CB, e’, always compete with holes in the VB, h". The overall 
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192 conductivity behavior depends on at the number densities of e’ and h‘, designated as n’ and n‘, and the 

193 mobilities, designated as p’ and p". If the product n’p’ is larger than n‘p‘, electrons control the overall 

194 conductivity and the semiconductor is said to be n-type. Conversely, if the product n’p’ is smaller 

195 than n'p‘, holes control the conductivity and the semiconductor is said to be p-type. 

196 n’p’>n'p‘ n-type [10a] 

197 n’p’<n'p‘ p-type [10b] 

198 Any semiconductor that is chemically pure is said to be intrinsic. In the intrinsic case, when electrons 

199 are thermally promoted across E g , the number of electrons on the CB will be equal to the number of 

200 holes in the VB: n’ = n*. Since the mobility of electrons on the CB is always higher than the mobility 

201 of holes in the VB, any intrinsic semiconductor will be n-type. While n’ and n* depend exponentially 

202 on the temperature, the temperature dependence of the mobilities, p’ and p* is weak. 

203 If a semiconductor is doped with aliovalent impurities, it can be made more n-type or p-type 

204 depending on the nature and concentrations of the dopants. Aliovalent dopants mean impurities with 

205 valencies other than those of the regular lattice constituents of the semiconductor material. However, 

206 aliovalency can also involve native lattice constituents of the semiconductor in a different valence 

207 state. This aspect will play an important role in the discussion below. 

208 The electrical conductivity of a material is the sum of the conduction of each type of charge carriers or 

209 each type of defects, all acting in parallel: 

210 a = 2aj = SniZipi [11] 

211 where n, is the concentration of the i h type of charge carrier, Zj is the charge, and p, the mobility. 

212 Minerals and rocks are “dirty” systems with complex chemistries and defects of many different kinds 

213 always present over a range of number densities or concentrations. Each type of impurity or defect 

214 contributes to o, though usually not more than two types will dominate in any given temperature 

215 window under a given set of thermodynamic conditions. 

216 

217 Activation Energies, Charge Carrier Concentrations and Mobilities 

2 1 8 Figures 2a/b show Arrhenius plots for mafic and ultramafic igneous rocks samples that have been 

2 1 9 conditioned by cycling them through the temperature interval of interest under controlled redox 

220 conditions, using a variety of heating and cooling rates [. Parkhomenko , 1967; Parkhomenko , 1982], 

221 From below 100°C to about 600-700°C, a for all rocks under study is consistent with eq. [8], i.e. with 

222 a single conductivity mechanism characterized by a single value of E act . Remarkably E atc is closely the 

223 same for all rocks, ~1 eV. This suggests that, irrespective of the type of rock and in spite of 

224 mineralogical and chemical differences, the conductivity of these rocks appears to be controlled by the 

225 same conductivity mechanism over the low temperature region (LT), from below 100°C to 600-700°C. 

226 Different samples of the same rock produce straight sections of the Arrhenius plots that can be shifted 

227 parallel to each other on the log a scale by a significant amount. Hence, in the AT window shown 

228 here, the only difference is the numerical value of a 0 . According to eq. [7] the only T-dependent 

229 contributions to a 0 are n and p, the number density of charge carriers and their mobilities, respectively. 

230 Since the temperature dependence of the mobility is weak, we can treat p’ and p" as approximately 

23 1 constant across the temperature interval of interest here, Hence our focus can center on n’ and n‘. 

232 
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Figure 2a/b: Electrical conductivity of mafic and ultramafic rocks after [Parkhomenko, 1982], 

Figure 2b shows that, when the rock samples were heated beyond 550-600°C, up to about 1000°C, 
straight sections with different slopes are obtained, suggesting different E ac t values for the high 
temperature (HT) region. 

The two linear sections points to two conduction mechanisms, endowed with activation energies, Ei in 
the LT range and E 2 , E 3 , E 4 . . . in the HT range, as suggested by eq. [9], At HT ionic conductivity sets 
in, exhibiting different slopes for different rock types, suggesting different ionic conduction with 
different E atc values. At the same time E atc in the LT range is closely the same for all rocks, although 
the 00 values may vary over orders of magnitude, even for the same rock type. For instance, the 
Arrhenius plot for peridotite in Figure 2a shows a conductivity over the 50-300°C window that is 6 
orders of magnitude higher than the conductivity of the peridotite shown in Figure 2b over the 100- 
500°C window. 


The constancy of the activation energy E act ~ 1 eV across different rock types in the LT range is as 
remarkable as the widely different oo values within one rock type. 


Nature of the Charge Carriers 

MgO, a II- VI material, is not a semiconductor. It has a very wide band gap of about 8 eV, MgO is the 
most ionic among all alkaline earth oxides with mostly Coulombic interaction between Mg 2+ and O 2 
[Pacchioni et al., 1993], As a main- group element with the electronic configuration Is 2 2s 2 , Mg has 
only one stable oxidation state, 2+. Therefore, pure MgO should be a model insulator. 

Large MgO single crystals can be melt-grown by the carbon arc fusion technique [Abraham et al., 
1971], Crystals are available at nominal 99.9% purity grades, designated as 3N, either undoped or 
doped with various levels with divalent 3d transition elements Mn, Fe, Co, or Ni. Purer MgO crystals, 
up to purity grade 99.99%, designated as 4N, are also available in sizes up to several cm 3 . 

If MgO single crystals are well-annealed, their electrical conductivity at room temperature can be 
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below 10' 16 [Q^cm' 1 ]. However, when heated, MgO single crystals display a striking dc conductivity 
behavior as illustrated in Figures 3a-d for a high purity 4N-MgO crystal measured under the most 
stringent clean conditions and with a guard electrode to remove contributions from surface 
conductivity [Kathrein and Freund, 1983], 

(a) 


900 800 700 600 500 400 300 2 SOX 



900 800 700 600 500 400 300 250X 




(d) 


0.8 


1.2 1.4 1.6 

1 000/T [K] 


1.8 



1 000/T [K] 


Figure 3. Sequence of dc conductivity measurements of a highest purity grade MgO single crystal 
under the most stringent clean conditions, (a) shows an anomalous increase in conductivity during the 
first heating, which is not observed as shown in (b) when the crystal is reheated immediately after 
cooling. However, if the crystal is annealed for several hours, the anomalous conductivity reappears 
as shown in (c). Panel ( d) shows that this anomalously high conductivity is due to the activation of 
charge carriers that conduct electricity with an activation energy ~1 eV [Kathrein and Freund, 1983]. 
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268 Upon first heating, starting around 350-400°C, the conductivity of the MgO crystal increases by about 

269 6 orders of magnitude as shown in Figure 3a. At the heating rate used in this experiment, 10° min' 1 , 

270 the conductivity reaches a maximum around 600-650°C and then decreases, only to increase again 

271 above 700°C. Around 850-900°C the conductivity joins an Arrhenius straight section that describes 

272 the reversible HT conductivity, which is endowed with an activation energy of 2.4 eV. At very slow 

273 cooling rates the HT mechanism can be followed downward to below 500°C [ Lempicki , 1953], At any 

274 cooling rate the HT mechanism is eventually replaced by an LT mechanism, endowed with an 

275 activation energy of ~1 eV. The transition from HT to LT occurs around 650-700°C. In the HT 

276 regime the electrical conductivity of MgO is controlled by Mg 2+ diffusion via an Mg 2+ vacancy 

277 hopping mechanism [Sempolinski and Kingery, 1980], possibly augmented by O 2 diffusion [ Wuensch 

278 et al., 1973], The nature of the LT 1 eV mechanism has long remained undetermined. 

279 If the MgO crystal is reheated immediately after cooling as in the run depicted in Figure 3b, the 

280 conductivity curve “2 nd up” closely follows the LT conductivity behavior observed during first cooling 

281 “1 st down” and also during second cooling “2 nd down”. If, however, the MgO crystal is not reheated 

282 immediately after cooling but allowed to anneal for several hours or overnight at temperatures below 

283 ~150°C, Figure 3c shows a conductivity pattern, “3 rd up”, similar to that of the first heating, “1 st up” 

284 in Figure 3a. This conductivity pattern is reproducible except that, during repetitive heating-cooling 

285 cycles, the highest conductivity achieved around 650°C slowly decreases with the number of cycles 

286 and with the holding time at the highest temperature. Another subtle but reproducible feature is that, 

287 during cooling, the transition from the HT section to LT section shifts to ever lower temperatures the 

288 more often the MgO crystal has been cycled through HT and the longer it has been kept at 

289 temperatures >900°C. 

290 A similar conductivity pattern has been observed in studies with natural samples, both single crystal 

291 minerals and rocks [ Constable and Duba, 1990; Duba and Constable, 1993; Roberts et al., 1999], In 

292 those cases the high conductivity recorded during first heating in the LT range has been blamed on 

293 adventitious surface contamination by carbonaceous films. Surface contamination can never be totally 

294 ruled out during conductivity measurements, even under very well controlled conditions. However, 

295 conductivity measurements can be performed in such a way that surface contamination as a reason for 

296 the anomalous electrical conductivity behavior can be ruled out [Kathrein and Freund, 1983], 

297 Figure 3d shows what happens when the heating of the 4-N MgO crystal is intermittently reversed 

298 three times within the range of the anomalously rapid conductivity increase, at 540°C, 555°C, and 

299 590°C. After each intermediate cooling, the MgO crystal was immediately reheated. For every cycle, 

300 Arrhenius straight sections are obtained, all with the same slope corresponding to an activation energy 

301 of ~1 eV. Each time, during reheating, the conductivity eventually rejoins the main branch, marked by 

302 solid red circles. The main branch is characterized by an E ac t value around 3.2 eV. 

303 Because this particular experiment was conducted at heating and cooling rates of 20°C min' 1 , the main 

304 branch in Figure 3d is shifted towards higher temperatures relative to the curves in Figures 3a-c, 

305 which were conducted at 10° min 1 . Except for this experimental effect, which due to the thermal 

306 inertia of the system, the conductivity patterns are very similar. 

307 A conductivity pattern as depicted in Figure 3d is not consistent with surface contamination. It points 

308 to two processes working in tandem. 

309 Process (i) involves the activation of defects in the MgO structure, which are electrically inactive but 

310 release mobile charge carriers as the temperature increases. These charge carriers must be electronic. 

3 1 1 The activation begins at a threshold temperature, ~350°C at the 10° min' 1 heating rate or ~450°C at the 

312 20° min' 1 heating rate, and has an apparent activation energy barrier (Ei + E 2 ) ~3.2 eV. Process (i) is 
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313 consistent with the number n of mobile charge carriers increasing exponentially as the MgO crystal is 

314 heated. Process (i) takes place only after the MgO crystal has been annealed for some time, operating 

315 between 350-600°C during the 10° min' 1 runs or 450-650°C during the 20° min' 1 run. Around 600- 

316 650°C process (i) comes to an end, suggesting that all available charge carriers have been activated. 

317 Above 700°C the conductivity decreases, indicating that charge carriers are being destroyed. This is 

3 1 8 consistent with the aforementioned decrease in maximum conductivity value obtained during 

3 1 9 consecutive heating-cooling cycles. 

320 Once the charge carriers are activated, they remain active for a sufficiently long time to immediately 

321 conduct successive cooling-reheating cycles. If the MgO crystal is allowed to anneal, the charge 

322 carriers return to an inactive state, from which they can be reactivated by reheating to 350-450°C. 

323 Process (ii) relates to how much current is carried by n mobile charge carriers, where n is the number 

324 of charge carriers activated through process (i). The charge carriers require an activation energy E 2 to 

325 jump from one site to the next. Both processes, (i) with the combined activation energy (Ei + E 2 ) ~ 

326 3.2 eV, and (ii) with the activation energy E 2 ~ 1 eV, follow Arrhenius relations. To describe the 

327 system under study we can use a relation very similar to eq. [9]: 

328 a = a 0 exp [-(Ei+E 2 )/kT] [12] 

329 The sum (Ei + E 2 ), ~3.2 eV, applies to the main heating branch in Figure 3d. E 2 ~ 1 eV applies to the 

330 cooling-reheating straight sections in Figures 3a-d. The value for Ei is ~2.2 eV, that for E 2 is ~1 eV. 

33 1 In summarizing this section we can state with confidence that highest purity-grade 4N-MgO single 

332 crystals exhibit a complex dc conductivity behavior in the temperature window between ambient and 

333 800-900°C. The conductivity increase by at least 6 orders of magnitude between 300°C and 600- 

334 700°C is clearly at odds with the presumed model insulator nature of MgO. It points to either the 

335 presence of impurities that are not included in the chemical purity grade or to the presence of defects, 

336 which had not been identified in prior work. 

337 There is one more important observation: MgO crystals doped with divalent 3d transition metal 

338 cations, Mn 2+ , Fe 2+ , Co 2+ or Ne 2+ at various, even high concentration levels exhibit the same electrical 

339 conductivity pattern as described here for 4N-grade MgO [. Freund et al., 1993; Kathrein and Freund , 

340 1983], This suggests that the anomalous conductivity behavior is not caused by transition metal 

341 impurities nor by charge carriers that would propagate via the metal cation sublattice. 

342 Furthermore we can state with confidence that the observed conductivity behavior of the MgO single 

343 crystals and their 1 eV conductivity mechanism cannot be due to surface contamination. MgO crystals 

344 that were temperature-cycled in pure 0 2 show the same conductivity pattern, ruling out that the 

345 increased conductivity would be caused by carbonaceous surface films [. Kathrein and Freund, 1983], 

346 

347 Valence Changes in the Oxygen Anion Sublattice 

348 MgO was the first oxide material for which evidence was obtained that the valence of oxygen anions is 

349 not fixed at 2-, as is widely assumed, but can change to the higher oxidation state 1-. When ultrahigh 

350 purity Mg(OH) 2 was thermally decomposed to MgO + H 2 0, forming a dilute MgO-H 2 0 solid 

351 solution, Mgi_5(OFl)260i_2&, it released copious amounts of H 2 , 10,000 ppm, even though the 

352 transition metal impurity level was <5 ppm, equivalent to a purity grade >99.999% [Martens et al., 

353 1976b], Since Mg cations are fixed in their 2+ valence state, the only possibility for H + to be reduced 

354 to H 2 was to concomitantly oxidize oxygen anions from valence 2- to valence 1-. 

355 Follow-on studies with MgO single crystals revealed that the redox conversion takes place in the MgO 
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356 matrix involving OH“ pairs associated with and chargewise compensated by Mg vacancies. During 

357 the redox conversion two hydroxyl protons are reduced to H, forming H 2 , while two hydroxyl oxygen 

358 are oxidized to O”, forming a peroxy anion, O 2 2 " [Batllo et al., 1991; Freund and Wengeler, 1982], 

359 Using the Kroger-Vinck point defect designation 1 we can write this reaction as: 

360 [OH*V Mg HO*]x <*> [o“(H 2 ) Mg ] x [13] 

361 The defect on the right hand side can release its H 2 molecule into an interstitial site: 

362 [o**(H2) Mg ] x ^ [o**V Mg ] x +(H 2 )i X [14] 

363 where the right side now describes an Mg 2+ vacancy chargewise compensated by a peroxy anion, 0 2 2- . 

364 If the MgO crystals are small [Martens et al., 1976a] or if the time for outdiffusion is long, the H 2 can 

365 escape from the solid matrix. Such loss of H 2 would change Mgi_ 5 (OH) 2 aOi _26 into Mgj_§0. If 6«1, 

366 cation-deficient Mgi _50 may also be written as MgO with excess of oxygen, MgO 1 + 5 . 

367 MgO is intrinsically diamagnetic. Peroxy anions are diamagnetic, because the two O” are tightly 

368 coupled and spin-paired. Upon heating, the spin pairing loosens up and eventually breaks, causing a 

369 transition from a diamagnetic to a paramagnetic state and dissociation. Magnetic susceptibility studies 

370 have provided information about this transition [Batllo et al., 1991], Peroxy spin decoupling and 

371 dissociation take place in two distinct steps: (i) the peroxy bond loosens up, (ii) full dissociation takes 

372 place. Designating as I the peroxy defect as I and as II the transition state, where spins associated with 

373 the two O begin to decouple, we have: 

374 [o”V Mg ] x » [ 0 -V Mg 0 -]x [15] 

375 I O II 

376 During stage II, the wave functions associated with the two O” delocalize over many neighboring 

377 oxygen anion positions, causing anomalies in fundamental physical properties such as in the thermal 

378 expansion and in the refractive index [Freund et al., 1994], However, at this decoupling stage, there is 

379 no significant increase in the electrical conductivity. 

380 During stage II, the two O dissociate. Dissociation involves an electron transfer from a neighboring 

381 O 2 donates into the [O* V Mg 0*] x defect. The neighboring O 2 ”, which had donated the electron, 

382 thereby turns into an O , which is no longer bound to the Mg 2+ vacancy. It becomes an unbound defect 

383 electron on the oxygen anion sublattice. The Mg 2+ vacancy retains one CT, forming a V center, a 

384 point defect widely studied by electron paramagnetic resonance spectroscopy [Henderson and Wertz, 

385 1977; Marfunin, 1979], Designating the dissociated state with III we have: 

386 [0-V Mg O-] x O [0*V M g]‘ + Cr [16] 

387 II O III 

388 Figure 4 outlines the temperature windows for the transitions between defects I, II, and III. It 

389 summarizes the insight gained from a wide range of studies aimed at understanding the electrical 

390 conductivity of nominally high purity MgO single crystals. These studies have positively identified 

391 peroxy defects in the MgO matrix as the cause for this anomalous conductivity behavior. 

391 

1 V stands for vacancy; subscripts identify the site (except for oxygen sites, where subscripts are omitted); 
superscript prime, dot, and x designate single negative, positive and neutral charges, respectively, double prime 
and double dot designate double negative and positive charges; subscript i means interstitial; square brackets 
outline the essential parts of any given point defect. 
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Unbound 0“ constitute defect electrons in the O 2- sublattice to be designated h' [ Griscom , 1990]. The 
If are free to roam - as electronic charge carriers - through the oxygen anion sub lattice. Therefore, 
when the peroxy defects break up, releasing h‘ charge carriers, the electrical conductivity increases 
sharply as demonstrated in Figure 3d, rising by 6-7 orders of magnitude over a 250°C wide interval, 
from about 450°C to about 700°C. 


Close to 700°C the release of h’ charge carriers comes to completion, marking the end of dissociation 
of defect II. Above 700°C the conductivity drops. The reason is that another process sets in, namely 
Mg 2+ diffusion, which becomes thermally activated above ~600°C [Lempicki, 1953], Mg 2+ diffusion 
is essential for non-stoichiometric Mgi_ 6 0 to return to equilibrium. It requires annihilation of the 
excess Mg vacancies at the crystal surface, a process that proceeds by way of cation vacancy 
hopping [ Conrad and Yang, 2010; Sempolinski and Kingery, 1980; Wuensch et al., 1991]. 


900 800 700 600 500 400 300 250°C 



1000/T [K] 


Figure 4: Interpretative representation of the 
anomalous electrical conductivity pattern of 
highest purity MgO single crystals, supposedly a 
model insulator. Three temperature windows are 
identified, below 250°C, between 250-450°C and 
450-700°C, each associated with different states of 
peroxy defects in the MgO matrix, designated as I, 
II and III. The dissociation of the defect III 
uniquely controls the electric conductivity of MgO 
through the release of positive hole charge 
carriers, h\ here shown as O’ states. 


Thermodynamic Equilibrium and Non-Equilibrium Conditions 

If O is to be introduced into MgO under thermodynamic equilibrium conditions, highly oxidizing 
conditions are required [Wuensch et al., 1991], No MgO crystal grown from the melt of a carbon arc 
fusion furnace under highly reducing conditions should contain oxygen in the valence state 1-. 
Therefore, the presence of peroxy defects appears to be in disagreement with thermodynamics. 

However, the reaction described by eq. [8], the loss of H 2 as described by eq. [9], and the presence of 
oxygen in the valence 1- are not in conflict with thermodynamics. Nor is the presence of excess Mg 2+ 
vacancies in the non-stoichiometric Mgi_50. The number of Mg 2+ vacancies in the MgO matrix is not 
a reflection of the 0 2 fugacity at the time of crystallization but a consequence (i) of the concentration 
of H 2 0 dissolved in the MgO matrix under thermodynamic equilibrium conditions and (ii) of the solid 
state reactions described by eqs. [8] and [9], which occurred at temperatures too low for equilibrium to 
be maintained. 

When nanocrystalline Mgi_50 produced by the outdiffusion of H 2 according to eq. [14] was heated 
above 600°C, atomic O started to evolve, indicating disproportionation [Martens et al., 1976b]: 

0 2 2 ' => O 2 " + O [17] 

Through this reaction Mgi_50 evolves toward the 1:1 stoichiometry that a material like MgO with 
Schottky defects is supposed to have at thermodynamic equilibrium [Alfe and Gillan, 2005]: 

Mgi_50 => MgO + 6/2 0 2 [18] 


12 


422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 

433 

434 

435 

436 

437 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 


Confirmation of Positive Holes by Dielectric Polarization 

Electrical conductivity measurements do not provide the sign of the majority charge carriers. If we 
want to learn more about the nature of the charge carriers, it is advisable to supplement conductivity 
measurements with corroborating techniques. 

One option is to conduct Hall effect measurements, widely used to study semiconductors [Clayhold 
and B. M. Kerns , 2008], which provide information about the sign of the majority charge carriers and 
their mobilities. Another option is to measure the dielectric polarization P, a fundamental property of 
materials, given by P = 8 0 (c-1)/4ji E extj where e 0 and 8 are the permittivity of vacuum and dielectric 
constant, respectively, and E ext is the externally applied electric field [. Freund et al., 1989]. 

Figure 5: Schematic 
representation of the dielectric 
polarization of insulators, 
giving rise to the effective 
dielectric constant £ e ff, a 
fundamental bulk property. 
Different contributions are 
identified. The first three on 
the right, at high frequencies, 
are characteristic of perfect 
crystals. The 4 th and 5 th arise 
from defects, which either 
increase the local polarizability 
or create mobile charges that 
can move as a space charge in 
the externally applied E field. 

As illustrated in Figure 5 the dielectric polarization P to tai is composed of several contributions in 
different frequency ranges. Five contributions are characteristic of the bulk, P e i, Pbond> Pdipoie, PiocaF 
and P S pace> pl us a sixth contribution, P SU rface> which arises from the surface [ Kittel , 1980], At high 
frequencies, in the visible (VIS) region, 1 0 1 5 — 1 0 1 3 Hz, only electrons can respond to the rapidly 
alternating E field. Their contribution P e i defines the refractive index n. In the infrared (IR) region, 
below 10 13 Hz, atoms linked by bonds or forming dipoles can follow the E ex t field, converting energy 
into atomic vibrations. The sum (P e i + Pb 0n d + Pdipole) defines the dielectric constant, 8, linked to the 
refractive index n by the Lorentz-Lorenz relation n = fz. In perfect dielectrics P to tai will not increase 
further with decreasing frequency. In practice defects are always present, some of which may have 
dipoles that can align in the E ex t field, possess a local polarizability that is different from that of the 
surrounding matrix. Their contributions are combined in Piocal- When mobile charge carriers appear, 
they can form space charges that respond to the external E ext field. Their contribution is designated as 
Pspace- In the context of this paper only Piocal and P S p ace . are of interest. 

The sixth contribution, P SU rface> arises uniquely from the surface. When a dielectric generates mobile 
charge carriers, these charge carriers repel each other in the bulk and are driven toward the surface. 

The sign of the surface charge is determined by the majority charge carriers, their number density n 
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and mobility | x. According to eqs. [5a/b] the surface charge is positive if n’p’ < n'p‘, indicating p-type, 
and negative, if nV>ny indicating n-type characteristics [King, 1984 #102], 

Figure 6 illustrates the principle of the measurement, the protocol and representative results [Freund 
et al., 1989], Figure 6a shows the sample placed between a small gold sphere and an annular gold 
electrode. Connecting the annular electrode to ground and applying a voltage to the small electrode, 
generates a divergent electric field, which produces a force F z acting on the sample. In practice the 
sample is suspended from a microbalance into a small, cup-shaped furnace, which contains both 
electrodes. Figure 6b shows a temperature-time program with temperature steps, AT, typically in 5° 
or 10°C increments/decrements. During each AT step, at constant T, a voltage-time program is 
applied as shown in Figure 6c with typical ±AV steps of ±5, ±10, ±20 or ±50 V. 


(d) 



Figure 6: Principle and execution of a dielectric polarization measurement in a divergent electric 
field gradient (see text). 

In the case of a dielectric with contributions from P e i, Pbond> Pdipoie, Piocab and P S p ace „ Fz is 
proportional to the square of the applied voltage leading to parabolas as in Figure 6d. In the case of a 
dielectric with mobile positive charges in the bulk, which make a P SU rface contribution, the parabolas 
are skewed by the addition of a linear term, proportional to the density of the surface charge, as 
illustrated in Figure 6e. 

Bulk and surface contributions can be separated by forming linear combinations F^ = ^ (F ± F ) 

1 — ± ± — 

and Fa = ^ (F - F ) , where F and F are the forces measured with ± and - bias, respectively. 

Zj Zj Z-i Z-i 

Fv depends quadratically on the bias voltage and is given by: 
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[ 19 ] 


= ~f {[Pel + P bond + P dipole + P local + ^spac J V ) ^ exA^ 
vol 

Fa depends linearly on the bias voltage and is a unique function of the surface charge density: 

^ =-/ [Psurface^E^dV 


[ 20 ] 


Figures 7a/b show Fv(T) and Fa(T) plots of an MgO single crystal, here a 3N-MgO crystal, 99.9% 

2 _|_ 

basic purity grade doped with 200 ppm Co . F^T) and Fa(T) graphs of undoped 4N-MgO or of 
moderately doped 3N-MgO crystals with divalent 3d-transition metal cations are practically 
indistinguishable. The data were obtained during stepwise heating to T max = 600°C in N 2 atmosphere 
following the protocol depicted in Figure 6b with 15 sec polarization time at each AY step during 
each AT step. The dotted line in Figure 7a indicates how F^T) changes, if the polarization time is 15 
min at each AV step. The 15 min polarization data were obtained in a separate run over this narrow T 
interval. Figure 7b indicate the evolution of F A (T). The black triangles show F A going negative, if 
heating is continued above 600°C. These data were obtained in a separate run. 


(a) 



Temperature [°C] 



Temperature [°C] , A 


Figure 7a/b: Bulk polarization TV and surface charge F A show the evolution of the effective dielectric 
constant and the thermal activation of positive charge carriers in MgO during heating to 600°C, 
which is consistent with the 2-step activation process as described by eqs. [15/16]. 

In Figure 7a the F^fT) data are fitted to the dielectric constant s of MgO, 9.6 at room temperature 
[Young and Frederikse, 1973]. At the limit of 0 Hz s begins to increase around 250°C at 15 sec 
polarization time or already around 200°C at 15 min polarization time, e reaches a plateau around 
350°C but exhibits a second stepwise increase beginning around 450°C. The Fa(T) data in Figure 7b 
indicate a slightly negative surface charge of the MgO up to about 350°C, followed by a pronounced 
positive surface charge, indicating the onset of p-type behavior. The positive surface charge reaches a 
maximum about 450°C, then decreases and turns negative above 600°C, indicating n-type behavior. 

According to eq. [14] F^T) is the sum of the five bulk contributions to Ptotal- Only the last two, Piocal 
and P S pace, depend on the presence of defects. The stepwise increase in Fv(T) indicates that these 
defects are thermally activated. During the first increase, there is no change in Fa(T), i.e. no activation 
of mobile charge carriers. This means that the first Fj(T) increase, which begins at 250°C or 200°C, is 
due to the thermal activation of defects, which only affect the local polarizability, Piocal. 
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495 This is fully consistent with the transition of the tightly bound, spin-coupled peroxy defect [q ** V Mg ] x , 

496 designated state I, to state II, the spin-decoupled state [O* V Mg 0*] x as described by eq. [10], It is 

497 consistent with the absence of any significant increase in the electrical conductivity of the MgO in the 

498 temperature range below approximately 350°C as evidenced in Figure 3a. 

499 By contrast, above 350°C, the electrical conductivity of the MgO starts to increase as shown in Figure 

500 3a, indicating the activation of mobile charge carriers. This is consistent with the increase in Fv(T), 

501 due to a contribution of P spa ce and simultaneous P su rface contribution as reflected in the change in Fa(T) 

502 in Figure 7b. This is consistent with the break-up of the peroxy bond as described by eq. [11], i.e. 

503 with the transition of the spin-decoupled state [O' V Mg 0'] x to the dissociated state [O' V Mg ]‘ + O’ and 

504 the activation of mobile positive hole charge carriers h\ 

505 Dielectric polarization measurements during cooling show that F^T) is fully reversible under the 

506 conditions of the experiments, i.e. under slow stepwise cooling over many hours. The high temperature 

507 portion of F^(T), the trend from the high positive values around T max ~500°C toward negative values 

508 at ~600°C, is also reversible. Below ~500°C, however, the surface charge lingers on, causing Fa(T) to 

509 remain positive down to ambient temperatures. At room temperature the surface charge dissipates 

510 over the course of hours with about the same time constant with which the electrical conductivity a(T) 

511 recovered during annealing as discussed in the context of Figures 3a-c. 

512 

513 Positive Hole Charge Carrier Activation in Feldspars 

514 MgO is not the only material that contains peroxy defects and activates positive hole charge carriers at 

515 distinct temperature steps upon heating. Silicates also contain peroxy defects, probably in the form of 

516 peroxy bridges, 03Si°° Si03 or 03Si °° A103, replacing 03Si ° v Si03 or OsSi^ AKf bonds 

517 respectively, or linking [Si0 4 ] groups. 

518 In the case of MgO the introduction of peroxy defects was linked to the dissolution of H 2 0 and 

5 1 9 formation of hydroxyl pairs associated with Mg 2+ vacancies, OH , followed by the redox conversion 

520 described by eqs. [13/14], Feldspars are tectosilicates the structures of which consist of [Si,A10 4 ] 

521 tetrahedra 3 -dimensionally connected by 03Si ° v Si03 and/or 03 Si°A 103 bonds. If H 2 0 is dissolved, it 

522 will most likely hydrolyze the bonds forming as primary product hydroxyl pairs: 

523 0 3 Si A Si0 3 + H 2 0 O 0 3 Si /OH HO/Si0 3 [21] 

524 The hydroxyl pair on the right hand side of eq. [21] could undergo the same type of redox conversion 

525 as OH” pairs in MgO, leading to” 

526 0 3 Si /OH Ho/Si0 3 O 0 3 Si /00v Si0 3 + H 2 [22] 

527 Feldspars indeed exhibit the same F^/T) and Fa(T) behavior when heated. This is exemplified in 

528 Figure 8a/b for F^fT) of a sanidine single crystal from Volkesfeld, Eifel volcanic district, Germany 

529 [Bertelmann et al., 1985], and for F^/T) of an anorthosite from Larvik, Norway [Wanvik, 2000]. The 

530 sanidine, Ko,96Nao,i4AlSi308, was cut from a very large (40 cm) single crystal extracted from a leucite 

53 1 phonolite tuff [Bank, 1967], The anorthosite sample consisted of an 8 mm diameter disk, 2 mm thick, 

532 cut of the coarsely polycrystalline Larvik anorthosite rock. The sanidine and anorthosite were 

533 transparent or translucent and smoky gray due to the presence of unbound positive holes, which 

534 localize at low temperatures at 03Si 0 AlC>3 or 03A1 0 A1C>3 [Speit and :ehmann, 1982], Both bleached 
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upon heating beyond 500°C to a light gray color. 


In both cases F^T) increases in 2 steps, at 200-250°C and 400-500°C for the sanidine and at 100- 
150°C and 400-500°C for the anorthosite. The F^(T) response, which is not shown, is somewhat 
complicated by the fact that the unbound positive holes anneal out during first heating. Therefore a 
stable Fa(T) response can be obtained only after careful annealing during second and subsequent 
heating cycles [unpubl. results]. 
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Figure 8a: FjfT) of a sanidine single crystal from Volkesfeld, Eifel volcanic field, brought up by an 
explosive eruption producing a leucite phonolite tuff [Bank, 1967]. The FyfT) response as measured 
at different voltages is shown. Figure 8b: FyfT) of a coarsely crystalline anorthosite from Larvik, 
Norway [Wanvik, 2000]. Both feldspar specimens exhibit the characteristic FyfT) behavior that 
points to the presence of peroxy defects and their 2-step activation, leading to the release of positive 
hole charge carriers. 

The similarity of the Fs(T) response of sanidine and anorthosite to the FvfT) response of MgO is 
striking. Since we know for certain that, in the case of MgO, the 2-step increase in F^T) is due to the 

2- step break-up of peroxy defects in the MgO matrix and activation of positive hole charge carriers, 
we take the results presented in Figure 8a/b as indication that these feldspars also contain peroxy 
defects. 

If this is so, the question arises whether there are other materials, compositionally similar to feldspars, 
which are produced under conditions that are not conducive to the introduction of “water” in form of 
solute 0 3 Si-OH. Normal plate glass, produced from the melt at relatively low temperature and low 
relative humidity levels, is such a material. Typical compositional ranges of clear glasses in mol-% 
are Si0 2 63-81, A1 2 0 2 0-2, MgO 0-6, CaO 7-14, Na 2 0 9-15, K 2 0 0-1.5 with most others below 0.5%. 

Figure 9a/b show the Fv(T) and Fa(T) responses of plate glass heated to 500°C. Fv(T) exhibits a 
increasing trend consistent with the increase in the overall polarizability as a function of temperature. 
At the glass transition temperature T g F^fT) increases more sharply, due to an increase in Pi 0C ai as a 
result of the softening of the glass structure and increase in diffusivity of Na + and other alkali cations. 
Fv(T) is mostly reversible under the conditions of the experiment with some hysteresis during cooling 
below T g . The reason for the slight deviation from a linear Fv(T) increase/decrease around 200°C is 
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unknown. However, the overall Fy(T) response of the glass is distinctly different from the F^T) 
behavior of either feldspar minerals shown in Figure 8a/b. There is no evidence for the characteristic 
F 2 (T) response linked to the 2-step activation of peroxy defects. 


(a) 



(b) 



Temperature [°C] 


Figure 9a/b: Dielectric polarization response of common plate glass provides no evidence for peroxy 
defects and the thermal activation of positive hole charge carriers. 

This conclusion is supported by the F^(T) response shown in Figure 9b. During heating a slight 
negative surface charge develops, which changes abruptly to a positive surface charge above T g . The 
evolution of F^(T) towards positive values is due to the high diffusivity of Na + and other alkali cations 
in the externally applied divergent E field. The Na + segregation is irreversible during cooling down to 
room temperature as demonstrated by the persistence of the positive F^(T) values in Figure 9b. 


The absence of evidence for peroxy defects in the dielectric polarization response of normal plate glass 
is consistent with the overall theme of this study: Peroxy defects are formed metastably out of solute 
OH”, CfSi-OH, or O3AI-OH, presumably hydroxyl pairs, which undergo a redox conversion during 
cooling. This reaction takes place in glasses as well as in crystalline materials. For instance, fused 
silica, which dissolves traces of H 2 O at a much higher temperature than ordinary plate glass, has been 
shown to contain peroxy defects and activated positive hole charge carriers in much the same way as 
melt-grown MgO [Freund and Masuda, 1991], 


Peroxy Defects and Positive Hole Charge Carrier Activation in Olivine 

The dielectric polarization behavior described in the preceding section indicates the thermal activation 
of mobile positive hole charge carriers in MgO and feldspar crystals but not in regular plate glass. It 
provides unambiguous support for the presence of peroxy defects in MgO grown from the melt under 
the extremely reducing conditions of a carbon-arc fusion furnace [Abraham et al., 1971] and in 
feldspar specimens that crystallized in a magma under reducing lower crustal conditions. The 
dielectric polarization response is fully consistent with the electrical conductivity behavior of the melt- 
grown MgO and magma-grown feldspar crystals. It is fully consistent with our conclusion from the 


18 


589 

590 

591 

592 

593 

594 

595 

596 

597 

598 

599 

600 
601 

602 

603 

604 

605 

606 


607 

608 

609 

610 
611 
612 

613 

614 

615 

616 

617 

618 


first part of this paper that positive hole are the controlling charge carriers that define the electrical 
conductivity in the LT region, up to about 600°C, of MgO and of mafic and ultramafic rocks. 

This study points to a mechanism by which crystals grown under reducing, even extremely reducing 
conditions, can acquire excess oxygen in the form of peroxy defects. The peroxy defects are 
introduced via the dissolution of H 2 O. The dissolution of H 2 O, however, is a process that is very 
common in nature. It takes place whenever minerals and rocks crystallize from an H20-laden magma. 

The electrical conductivity of dunite, an olivine-rich peridotite, is illustrated in Figure 10, adapted 
from work by Constable and Duba [ Constable and Duba, 1990], During heating in a 1:1 CO:C02 
atmosphere, the conductivity increases rapidly between ~350°C and 670°C. At 670°C, as T was held 
constant for an unspecified length of time, the conductivity decreased precipitously. During further 
heating the conductivity increased again, merging into a high temperature branch similar to the HT 
branch observed with MgO. During cooling the beginning of an Arrhenius straight section can be seen 
with an activation energy of ~1 eV characteristic of the LT conductivity response of MgO. 

The interpretation given by Constable and Duba [ Constable and Duba, 1990] is based on the 
assumption that, during heating, carbon was depositing out of the vapor phase onto the sample surface, 
providing a low resistivity pathway. The carbon deposition is assumed to have ended at 670°C when, 
according to thermodynamic calculations, the 1:1 CO:C 02 mixture was leaving the stability field of 
graphite, creating oxidizing conditions, which caused the surface carbon film to “bum off’. 


Figure 10: Electrical conductivity measured 
on a sample of Jackson County peridotite, 99% 
olivine, in a 1:1 CO: CO 2 atmosphere showing 
an anomalously high conductivity in the LT 
region, thought to be due to vapor-phase 
deposited surface carbon (after [Duba and 
Constable, 1993] f 


727 560 441 352 282 227 181 143°C 



The comparison between electrical conductivity and dielectric polarization shows that P to tai and its 
separation into Pbuik and P su rface can provide unique information about the nature of the charge carriers 
responsible for the conductivity response. We therefore present next the dielectric polarization 
response of an olivine single crystal. To circumvent complications arising from carbon precipitation in 
the graphite stability field we chose to conduct the first set of experiments in high purity N 2 . The 
choice of an inert gas atmosphere is justified on the basis of the observation that the CO:CC >2 =1:1 
mixture remained unequilibrated up to at least 500°C, possibly higher [ Constable and Duba, 1990], 

The sample was prepared by cutting a 1 mm thick slice, 8 mm diameter, from a randomly oriented 
cylinder of a gem-quality San Carlos olivine single crystal, 7 mm thick, polished on both sides. Before 
the slice was cut off, the IR spectrum was recorded. 

As shown in Figure 11 the broad absorption in the spectral range of the voh stretching vibrations, 
3150-3700 cm' 1 , lacks distinct features except for a weak band at 3600 cm -1 . The absence of voh 
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bands is conventionally interpreted to mean that such an olivine crystal without OH” has a very low 
solute “water” content and may have come from a “dry” of the upper mantle [Khan and Shankland, 
2011a], However, there is also the possibility that, similar to OH” pairs in MgO, solute OH” in olivine 
may have converted to H 2 plus peroxy defects. In this case the near-absence of voh bands may mean 
that an olivine crystal with an IR spectrum as exemplified in Figure 11 may have lost most of its 
solute OH” content due to this redox conversion. If the solute OH” content was originally high, the 
olivine should have a high concentration of peroxy defects. Measuring the dielectric polarization 
would then be the technique of choice [Freund et al., 1989], 


Figure 11: Infrared spectrum of a 7 mm thick 
section of a gem-quality olivine single crystal in 
4000 3800 3600 3400 3200 3000 the range of the O-H stretching frequencies 

Wavenumbers [cm ] indicates an extremely low solute OIF content. 

Figures 12a/b show Fv(T) and Fa(T) of this olivine crystal, recorded according to the protocol 
outlined in Figure 5 with 15 sec polarization time at each AT and AV step. F^T) is fitted at room 
temperature to a dielectric constant 8 scale. Because the orientation of the olivine single crystal disk is 
not known, we chose for the one-point fit the average of the three components of the dielectric 
constant along the crystallographic axes a, b, and c of forsterite [Cygan and Las aga, 1986], 

(a) (b) 





Temperature [°C] Temperature [°C] 


Figure 12a/b: Bulk polarization Fy and surface charge F ^ parameters of a San Carlos olivine single 
crystal, which show that, despite the chemical complexity of this mineral, its effective dielectric 
constant increases upon heating in the same way as that of high purity MgO, in two distinct steps, and 
a positive charge appears around 400 °C. 

Fv(T) of this olivine crystal is strikingly similar to Fv(T) of MgO and of the two feldspar specimens. 
During heating the bulk polarization increases in a stepwise fashion at 100-250°C and 450-500°C. 
During cooling F^T) follows closely the heating curve down to at least 100°C. 
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639 Fa(T) recorded is also strikingly similar to Fa(T) of MgO and the feldspars. The curves labeled +/- 

640 and -/+ correspond to Fa values obtained during heating taking the data from the sections going in AV 

641 steps from +500 V to -500 V and from -500 V to +500 V, respectively. The differences between the 

642 two F a (T) curves provides information about the time-dependent polarization during the 15 sec time 

643 steps. Fa(T) remains close to zero from 100°C to almost 400°C, while Fy(T) goes through its first 

644 stepwise increase. Thus the first F^T) increase is due to a Pi 0C ai contribution to Pbuik, without 

645 activation of mobile charge carriers. Mobile charge carriers only appear upon approaching 400°C. At 

646 first their number is small, insufficient to cause a noticeable increase in F^T). By about 450-500°C 

647 their number has increases enough so that F^T) also increases, due to a P sp ace contribution to Pbuik- 

648 The sign of the mobile charges in the olivine is unambiguously positive, as in MgO and the feldspar 

649 specimens, providing firm evidence that positive holes are the active charge carriers. 

650 

65 1 Metastability and Thermodynamic Equilibrium 

652 The very weak IR absorption of a 7 mm thick olivine crystal in the range of voh stretching vibrations 

653 as evidenced by Figure 11 indicate a near-complete absence of solute OFT. The same olivine crystal 

654 however provides evidence for positive hole charge carriers that become thermally activated at about 

655 the same temperature at which positive holes are activated in MgO, suggesting that they derive from 

656 peroxy defects in the olivine matrix, formed by the same type of redox conversion that turned OH in 

657 MgO into O 2 2 + H 2 as given by eq. [13]. For olivine we have: 

658 0 3 Si/ OH #HO/Si0 3 O 0 3 Si/°°\Si0 3 + H 2 [23] 

2 _|_ 2 1 

659 where the # sign is used to indicate a divalent cation vacancy, Mg or Fe . 

660 By analogy to the dilute solid solution MgO-H 2 0, we formulate a dilute olivine-H 2 0 solid solution as: 

661 (Mg,Fe 2+ ) 2 Si0 4 + 6 H 2 0 O (Mg,Fe 2 + ) 2 - 6 Si 04 - 2 &OH 2 6 [24] 

662 If the OH“ in this olivine-H20 solid solution undergo a redox conversion according to eq. [23] and if 

663 some or all of the molecular H 2 is lost from the crystal matrix by outdiffusion, we have: 

664 (Mg,Fe 2+ )2-&Si0 4 -2&0H 2 s O (Mg,Fe 2+ ) 2 _6Si0 4 + 6 H 2 [25] 

665 where we do not separately indicate on the right hand side the fraction of oxygen in the valence 1-. 

666 We specify, however, that the iron continues to be ferrous, Fe 2+ . In order for the cation-deficient 

667 olivine on the right hand side to be balanced, a 26 fraction of the oxygens have to be in the valence 1-. 

668 Low solute OH“ concentrations are usually interpreted to mean that such crystals come from “dry” 

669 regions of the upper mantle [Bell and Rossman, 1992; Miller et al., 1987; Rossman, 1996], Eq. [25] 

670 suggests instead that solute OH may have “disappeared” from the IR spectra by converting to H 2 plus 

671 peroxy. Molecular H 2 , if present, has an extremely weak vh-h absorption [ Kriegler and Welch, 1968]. 

672 Other techniques such as SIMS or nuclear reaction profiling, which measure the total H concentration, 

673 often yield higher H values than those obtained by IR measurements of voh [Bell et al., 2003]. 

674 

675 “Burning off’ Positive Holes in Olivine 

676 if the enhanced electrical conductivity observed during initial heating were due to carbonaceous films 

677 precipitated within the stability field of graphite [Constable and Duba, 1990] ? such films should “bum 
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off’ under oxidizing conditions. If the enhanced electrical conductivity is due to h‘, i.e. highly 
oxidizing O in a matrix of O , the best way to annihilate them is to use reducing conditions. 

We therefore conducted a dielectric polarization experiment with a San Carlos olivine single crystal, 
heating it first in pure N 2 atmosphere in 20 temperature steps from 50°C to 600°C, cooled it back in 
the same steps, and let is sit overnight in the pure N 2 . F^CO as recorded during this first heating is 
shown by the top curve labeled with squares in Figure 13a. 

Next we replaced the N 2 with a 9:1 N 2 -CH 4 gas mixture and reheated the same olivine crystal five 
times. The 1 st heating in N 2 -CH 4 , up to 575°C, produced the F 2 (T) curve marked by open circles in 
Figure 13a. After stepwise cooling the olivine crystal was reheated in the same N 2 -CH 4 gas mixture 
four more times. The 3 rd heating up to 550°C produced the Fy(T) curve shown in Figure 13a by open 
triangles. The 5 th heating up to 600°C produced the Fy(T) curve in Figure 13a marked by diamonds. 

Figure 13b shows Fa(T) curves recorded during the 1 st heating in the N 2 -CH 4 gas mixture with ± bias 
voltages of 20V, 40V, 60V and 80V. The Fa(T) response displays the same overall pattern as observed 
for MgO and the feldspar specimens. The dashed lines present the Fa(T) response during stepwise 
cooling: fully reversible from 600°C to about 500°C, followed by elevated Fa(T) values below 400°C 
due to the fact that, after heating in the N 2 -CH 4 mixture, the positive surface charge no longer 
dissipates after cooling. 

Comparing the F 2 (T) response during the 1 st run in pure N 2 and the subsequent runs in the N 2 -CH 4 gas 
mixture provides insight into the reactions at the olivine surface. Figure 13a shows that, during 
heating in the presence of CH 4 , the Fv(T) response i s similar up to about 450°C, though the values are 
slightly lower. This overall decrease of F 2 (T) is due to the fact that the 1 st heating in pure N 2 had 
already caused changes in the polarization response of the olivine, presumably due to oxidation of Fe 2+ 
to Fe 3+ in the olivine matrix, equivalent to the annihilation of mobile h* charge carriers. 


701 Above 450°C the F 2 (T) curve drops during heating in the N 2 -CH 4 gas mixture, indicating that h' 

202 charge carriers are being consumed. Above 530°C Fz(T) increases again, indicating that the surface 

703 reaction, which consumes h‘ charge carriers, has stopped and new h‘ are able to flood the surface. This 

704 behavior, in particular the dip of F 2 (T) above 450°C, repeats itself during the 3 rd heating cycle. By the 

705 time of the 5 th heating, however, the dip has essentially disappeared, indicating that the reaction 

706 between the olivine and CH 4 has reached some form of endpoint. 


707 Of further interest is the fact that, in the AT range marked in pale yellow in Figures 13a/b, Fv(T) 

708 continues to increase even though the Fa(T) curves show that h' charge carriers have already begun to 

709 flood the olivine surface. Fv(T) begins to dip in the AT range marked in yellow, i.e. only after the 

710 F A (T) values have reached their maxima. The area marked in red outlines the AT range, where F 2 (T) 

711 recovers, rising in spite of the presence of CH 4 . 
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Figure 13a: Bulk dielectric polarization of an 
olivine crystal subjected to 6 heating /cooling 
cycles, the 1 st in pure N2, the others in N2-CH4 
gas mixtures. The colors indicate temperature 
regions at which K become activated but do 
not yet react with the CH4 (light yellow), start 
to react with the CH 4 (yellow), and overwhelm 
the reaction with CH4 (red) 


Figure 13b: Surface charge developing on an 
olivine crystal during the 1 st heating/cooling in the 
N2-CH4 gas mixture. In contrast to the first heating 
in pure N2, when the positive surface dissipates after 
cooling, heating in the N2-CH4 gas mixture leaves a 
positive surface charge. The colors correlate with 
AT ranges at which the bulk polarization of the 
olivine crystal is affected by the reaction with CH4. 


Summarizing this last set of experiments we can say with additional confidence that, though the 
olivine crystals under study have originated in the upper mantle and crystallized in a highly reducing 
environment, they release highly oxidizing h‘ charge carriers upon heating. From this we infer that 
they contain peroxy defects, presumably formed during cooling by way of a redox conversion of OH 
pairs that had been incorporated into the olivine structure at upper mantle temperature through 
dissolution of traces of “water”. Consistent with our earlier conclusion that h* are highly oxidizing O , 
they are consumed by a reaction with CH4, presumably through oxidative coupling, which produces 
OH and methyl as an intermediate, *CH 3 [Asame et al. , 1988], 


Discussion 

In this paper we have focused on fundamental processes that are important for understanding the 
electrical properties of materials, both laboratory-grown and from natural environments, both single 
crystal minerals and igneous rocks. We make the extraordinary claim that MgO crystals, grown from 
the melt under the viciously reducing conditions of a carbon arc fusion furnace, contain peroxy defects 
in their crystal structure, hence oxygen in the oxidized valence state 1-. When the peroxy defects 
break up, they release positive hole charge carriers, formally defect electron in the oxygen anion 
sublattice, equivalent to O in a matrix of O 2 . These positive holes have two outstanding properties: 
they are highly mobile and highly oxidizing. 

We make the even more extraordinary additional claim that olivine crystals, which have been brought 
up from the upper mantle, feldspar crystals, which have been brought up from the lower crustal depth, 
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732 and - more generally - mafic/ultramafic rocks such as basalts, gabbro and peridotite produce the same 

733 type of mobile positive hole charge carriers and, hence, must also contain peroxy defects. This claim 

734 seems to run counter to the well-established fact that upper mantle and lower crustal environments are 

735 highly reduced and highly reducing. Thermodynamically, oxygen in the oxidized valence state 1- can 

736 never be introduced under these in situ conditions. 

737 In fact there is no contradiction with thermodynamics. The reason is that the process, which occurs 

738 under in situ conditions, is the dissolution of H 2 0 into the matrix of MgO grown from the carbon arc 

739 fusion melt and into the matrix of silicates crystallizing from a magma. The dissolution of H 2 0 and 

740 formation of solute OH or 0 3 Si-0H are redox insensitive. They can occur irrespective of the redox 

741 state of the system. According to our results O is introduced through a follow-on process, the redox 

742 conversion of hydroxyl pairs to peroxy defects plus H 2 . This process only occurs under metastable 

743 conditions, namely during cooling when the systems have drifted out of thermodynamic equilibrium. 

744 Depending upon the rate at which a rock or a crystal is cooled, thermodynamic equilibrium can only 

745 be maintained as long as diffusional processes remain thermally activated. Diffusional processes are 

746 necessary for the adjustment of the oxidations states of transition metal cations and the adjustment of 

747 the cation vacancy concentration, both dictated by the oxygen fugacity of the system. Without these 

748 continuous adjustments any system will be out of thermodynamic equilibrium. Even at very slow 

749 cooling rates these diffusional processes become so slow around 600°C that they can be considered 

750 frozen 2 . At faster cooling rates the freezing temperatures tend to be higher, more in the 700°C range. 

751 The in situ redox conversion, which turns solute hydroxyl pairs into peroxy defects plus molecular H 2 , 

752 takes place below 600°C, in the 600-400°C window, when the systems under consideration - MgO, 

753 feldspar, olivine crystals and igneous rocks - have already drifted out of thermodynamic equilibrium. 

754 Since the redox conversion as outlined by eqs. [13] and [23/25] operate outside equilibrium, oxygen in 

755 the valence 1- can exist metastably in an otherwise reduced environment and can co-exist with 

756 reduced cations such as Fe 2+ . Since O form highly mobile electronic charge carriers, the electrical 

757 conductivity behavior is affected. 

758 Therefore, the claim of oxygen in the valence state 1- is not as extraordinary as it may first appear. 

759 The presence of O does not contradict thermodynamics because, by the time O is introduced, 

760 equilibrium thermodynamics is no longer in control. The work presented here casts doubts on the 

761 validity of the widely offered explanation of thin carbonaceous films deposited from C0/C0 2 -laden 

762 fluids in the crust [Glover and Vine, 1992; Shankland et al., 1997] and upper mantle [Constable and 

163 Duba, 1990; Duba and Constable, 1993; Mathez, 1987; Roberts et al., 1999], Maybe electrolytically 

764 conductive, intergranular fluids and/or intergranular carbonaceous films are not needed to explain the 

765 origin of zones of high electrical conductivity in the lower continental crust, which have been called 

766 “a long-standing mystery” [Frost et al., 1989], Those high conductivity zones, which have been 

767 identified by magneto-telluric sounding techniques, appear to be an almost universal feature of the 

768 Earth’s crust at a depth range of about 20 km [Glover, 1996; Lee et al., 1983; Qin et al., 1992], 

769 If peroxy defects and thermally activated positive hole charge carriers can be metastably present, we 

770 can tie the high conductivity zones to the geotherm [Freund, 2003], If temperatures along the 

771 geotherm reach values sufficiently high to thermally dissociate peroxy defects and release positive 

772 hole charge carriers, we can envision a depth range, where the peroxy defects would be in a 

773 permanently dissociated state. With reference to eq. [7] this means an increase in the number density 

774 n of mobile charge carriers and, hence, an increase in the conductivity a. Deeper into the crust, at 

774 

2 

This statement holds for “dry” systems. In environments with free supercritical fluids and other fluid- 
dominated systems equilibrium can be established down to lower temperatures. 
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775 higher temperatures, the positive holes will be consumed by equilibrium oxidation reactions, for 

776 instance Fe 2+ to Fe 3+ , concomitant with structural rearrangements that require cation and cation 

777 vacancy diffusion to be thermally activated. 

778 An analogous situation is seen in the experimental conductivity curves obtained during heating of 

779 MgO single crystals as depicted in Figures 3a, 3c, and 3d: at temperatures above 500-700°C the 

780 electrical conductivity decreased precipitously because the h* charge carriers begin to be consumed by 

781 reactions linked to the return of the MgO crystals toward equilibrium. The same behavior is reflected 

782 in the conductivity curve recorded during heating of an olivine single crystal as depicted in Figure 9, 

783 even though the authors of this graph base their interpretation of the conductivity decrease at 630°C on 

784 the assumption that vapor-phase deposited carbon films on the sample surface where “burning away” 

785 as the system moved out of the stability field of graphite in the C0:C0 2 gas mixture [ Constable and 

786 Duba, 1990], Finally, the dielectric polarization behavior of the olivine single crystal heated in a 

787 CO:CH 4 gas mixture as depicted in Figure 13a/b indicates as well that h* charge carriers are being 

788 consumed under highly reducing conditions - in this case by reacting with CH 4 by way of an oxidative 

789 coupling that produces OH and methyl as an intermediate, *CH 3 [Asame et al., 1988], 

790 
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